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Abstract 

We prove some regularity results for viscosity solutions to strongly de- 
generate parabolic semilinear problems. These results apply to a specific 
model used for pricing Mortgage-Backed Securities and allow a complete 
justification of use of the classical Ito's formula. 
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1 Introduction 

In this paper we investigate the second order regularity of bounded viscosity 
solutions to the following semilinear parabolic equation of degenerate type 

d t u + H{x,t,u, X7u, V 2 u) = 0, (1.1) 

where the hamiltonian function is defined by the expression, 

H(x,t,u,p,X) = -^tr{a* T (t)X) + (n(x,t),p)+\(u)\* T p\ 2 

+ V (u)(a T (t)p,w(x,t)) + f(x,t,u), (1.2) 

for every (x,t) e M. N x (0, T), a < u < b, p E R N , X e S N , where S N is 
the space of iV x N symmetric matrices endowed with the usual ordering. In 
(1.2), tr, | • | and (,) denote the trace of a square matrix, the Euclidean norm 
and inner product, respectively. Moreover if AiNxdi^-) is the space of real 
N x d matrices, with N > d, then we assume that, /i:l"x [0,T) -> R N , 
a : [0,T) -> MwxdW, w:R N x [0,T) -> R d , and / : x [0,T) x (a, b) -> R, 
are continuous functions. 

Actually our main motivation comes from the following semilinear equation 
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d t U-^tr(aa T V 2 U)-(^U)+ PT ^^ } +r(U + h)-Th = 0, (1.3) 

in R N x [0, T), p > 0, t, T > 0, where £ = and /i = i) are regular func- 
tions of their variables. Equation (1.3) has been proposed in [31] as a differential 
model used for pricing some widely traded American financial instruments, the 
Mortgage-Backed Securities (MBS); in particular the model was derived fol- 
lowing the outline of X. Gabaix in [13]. 

Let U be a viscosity solution of (1.3), then setting u = U + h + u solves a 
differential equation of type (1.1). The arbitrage pricing principle applies to 
financial instruments whose cash flows are related to the values of some eco- 
nomic factors, such as the interest rates (r = r(t)). Using that principle the 
value a these securities can be expressed as a conditional expectation over the 
probability space (f2) of the underlying factors that determine the instrument's 
price (Vt) with respect to a particular measure (Q) defined over that space. As 
a consequence the knowledge of second order regularity properties of u, and 
therefore of U, plays a fundamental role in order to close in a rigorous way the 
financial argument used for deriving equation (1.3). If we have sufficient condi- 
tions about U, for applying the classical Ito's rule (see for istance [23], or [1]), 
then, as we will show in section 5, there exists a probability measure called the 
Equivalent Martingale Measure, Q, which has the required empiric and statistic 
representation proposed in [13], such that the following probabilistic equality 
holds: 

U{X t ,T-t) = E t[J t (T-r(T-s))e-f° r ( T - K)dK h{X s ,T-s)ds\, a.s.(1.4) 



dX t = n(X t ,T-t)dt + cr(T-t)dWt, < t < T. (1.5) 

There ¥,f [■] denotes the conditional mean value up to the time t, taken with 
respect to Q, and £(t) = eti r ^ ds , U(-,0) = 0, h > 0, h{-,0) = 0, r is a deter- 
ministic function which represents the interest rate variable. The process Wt is 
a d-dimcnsional standard Brownian Motion over the probability space fi, while 
X t is a A-dimensional Ito process which describes the factors which affect the 
value of a Mortgage-Backed security. In this model h models the remaining 
principal during the life of the pool of mortgages (see [13], [33]). 
The equation (1.4) represents the conclusive statement of our works about the 
pricing equation introduced in [31] and then studied from the point of view of 
the existence and uniqueness. Actually that formula gives the formalization of 
the existence of the market price of risk proposed by X. Gabaix in [13], which 
describes the martingale measure such that the value of the security can be 
expressed as the solution of the Hamilton Jacobi equation (1.3). 
Our main result concerns the semiconvexity /concavity property, at a fixed time 
t G [0, T), of the viscosity solution u of the equation (1.1). Applying our compar- 
ison result presented in [32] , it can be proved the existence of a unique viscosity 
solution for the equation (1.1), with an assigned continuous and bounded initial 
datum uo, valued in the interval (a, b). Moreover, using the same arguments of 
Theorem 4.5 in [32], it can be also deduced that, if the coefficients /x, A, 77, / 
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and the initial datum are i-uniformly Lipschitz continuous functions then, also 
the viscosity solution has the same regularity at a fixed time. Keeping in mind 
these features, we start our considerations from the assumption of the existence 
of a t-uniformly Lipschitz continuous solution u for the equation (1.1) and, then, 
we prove the semiconvexity /concavity property of it. 

The first result shows that u(t) is a semiconvex function over the whole space. 
In the recent literature of viscosity solutions for Hamilton- Jacobi equations, this 
kind of property is not still proved for a class of problems which could include 
our differential model. In [21], the authors only prove the semiconcavity of 
the solution for the Bellman equations, but the type of nonlinearity of their 
equation, does not include our Hamiltonian (1.2). There is another important 
paper, by Y. Giga, S. Goto, H. Ishii, and M.H. Sato, [14], where it was proved 
the convexity preserving property of the solution. But also in that case the au- 
thors study a fully nonlinear equation, whose Hamiltonian does not depend on 
the unknown u. Therefore, although the initial datum in the original financial 
model (1.3) is constant, their technique is not suitable for treating our equations 
(1.1), (1.3). 

Using our results about the semiconvexity and semiconcavity and assuming the 
same regularity of the coefficient used in (2), we deduce the W 2 '°° regularity, 
uniformly in t G [0, T). 

In Theorem 2.5, we conclude next, with a global result of regularity which con- 
tains also the time regularity of the solution. In particular, this last result will 
follow through an application of the Lipschitz continuity of u with respect to 
the time variable, which we will present in section 4. 

The general framework for deriving a financial princing equation consists in the 
application of the classical Ito's rule. Therefore, we should consider classical so- 
lutions which have continuous derivatives. However, there are some works about 
some generalizations of that formula, see for istance in [17]. Here it is proved 
that the Ito's formula holds in arbitrary dimensions for / € W 2 ' 1, °°(IR Ar x (0, T)), 
if the equation 

df(X t ,t) = Cf(X t ,t)dt + V T f(X t ,t)a(t)-dW t , a.s., (1.6) 

where 

Cf(x,t) = d t f(x,t) + (n(x,T-t),Vf(x 1 t)} 
+ ±tr(aa T (T-t)X7 2 f(x,T-t)) 7 

are interpreted appropriately using the generalized Hessian. 
In section 5 we derive some mathematical and financial consequences concerning 
the model (1.3). In particular we present a result which combines the regularity 
of the solution of the model (1.3) and the results of [17] to state the relation 
(1.4). To reach this objective we propose a technical condition about the pro- 
cess X t which allows us to extend the Ito's rule as a conditional expected value 
representation, see Lemma 5.4 later. The technical results proposed in section 5 
allow to consider the degeneration of the process X t , produced by the inequality 
N > d. That degeneration was studied in literature to give sufficient conditions 
to state the existence of densities for the solutions of stochastic differential equa- 
tions, see for instance [5], which contain a probabilistic form of Hormander's. 
Our condition is of a different nature and applies to many practical cases, like 
for constant coefficients, where the resultsa of [5] do not hold. 
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2 Main Results 



In this section we present our main results about the regularity of the solutions, 
the Section 3 being devoted to their proofs. 

Here wc leave out the definition of a viscosity sub/super solution of a differen- 
tial problem like (1.1), and refer the reader to some classical works about the 
viscosity theory, such as [10] or [25]. 



Definition 2.1. A function g in R w is said to be semiconvex with constant 
L > 0, if 

g(x + h) + g(x - h) - 2g(x) > -L\h\ 2 , (2.7) 

for every x, h € M. N . 

In a same way we will say that g is semiconcave with constant L > if — g 
is semiconvex with constant L. 

Over R m , we will denote (•, •) the usual standard scalar product, while, over the 
space, A4 m , n (R), of matrices with real coefficients and m-rows and n-columns, 
we will consider the usual norm, 

||A|| = sup \Ax\, VAeM m . n {R). 

Moreover we will denote as Im(A) the range of the linear map defined by the 
matrix A. 

If g : M. N x [0, T) x / — > R, where / is a closed interval in R, then we will use 
the following notations, 

IMloo = esssup( x ^ u ) eR N x[0iT ) XI \g(x,t,u)\, 



\\g(t)\\oc = esssup (xu)eR N XI \g(x,t,u)\, Vie[0,T). 

We will denote as W fc <°°(IR w ), k = 1,2, the usual Sobolcv Space of bounded 
functions with weakly bounded derivatives in R N of order less o equal to k. 
Moreover W 2 ' 1 '°°(IR Ar x (0,T)) denotes the Sobolev space of functions u, with 
weakly derivatives d t u, d l u,df J u e ^^(R^ x (0,T)), for i,j = 1, . .. ,N. These 
spaces are respectively endowed with the following norms 

k 

i=0 

IMIws.loo = HftuHoo + ||V«||oo + ||V 2 u||oo. 

If the functions take values in R M , we refer the same properties to the single 
components, and a similar notation holds for functions which depend on another 
variable in /. 



In the following we shall assume the existence of a viscosity solution u of 
(1.1), such that u(-,t) e W^°°(R N ) and whose norm is uniformly bounded in 

te [o,T). 
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Theorem 2.2. Let the viscosity solution u of the equation (1.1) be valued in 
a bounded closed subinterval I of the domain (a, b), and let [J,(-,t), w(-,t) G 
W 2 >°°(R JV ) ; /(-,*,-) e W 2 >°°{R N x J), uniformly in time. 
Assume: 

i) \eC((a,b)), V eC((a,b))nC 2 (I); 

ii) For every (x,t) G WL N x [0,T), io(ar,t) G Im(a T (t)). 

If uq is semiconvex, then there are positive constants C, M , Co, such that 

u(x + h,t) + u(x - h,t) - 2u{x,t) > -{M e ct + C )\h\ 2 , (2.8) 

holds for every x, h G R N , and t G [0, T). Therefore, for every t G [0, T), u(-,t) 
is semiconvex. 

We have also the analogous result for the semiconcavity of u. 

Theorem 2.3. Let the viscosity solution u of the equation (1.1) be valued in 
a bounded closed subinterval I of the domain (a,b), and let n(-,t), w(-,t) € 
W 2 '°°(R N ), /(-,*,-) € W 2 '°°(M W x I), uniformly in time. Assume: 

i) \eC((a,b)), V eC((a,b))nC 2 (I); 

ii) For every (x,t) G R N x [0,T), w(x,t) G Im(a T (t)). 

If uq is semiconcave, then there are positive constants C, M , Co, such that 

u(x + h,t) + u(x -h,t)- 2u(x, t) < {M e ct + Co)\h\ 2 , (2.9) 

holds for every x, h G R N , and t G [0, T). Therefore, for every t G [0, T), u(-,t) 
is semiconcave. 

As a consequence, we will deduce the announced result about the regularity 
of the solution u. 

Theorem 2.4. If fj,, w, f, A and r\ satisfy the assumptions of Theorems 2.2 
and 2.3 and u G W 2 ^{R N ), then u(t) G W 2 '°°(R JV ) ; uniformly in time. 

Theorem 2.5. If a is Lipschitz continuous, fi, w e W 2 ' 1 '°°(R JV x (0,T)), / G 

W 2,l,2,oo( R JV x (- 0jT ) x jj^ g / m ( CT T^) j or ewer2/ ^ ;f ) g R W x [Q )T ^ 

?i G W 2 '°°(R Ar ), i/ien, 

iie W 2 ' 1 '°°(R Ar x (0,T)). (2.10) 

To prove Theorems 2.2, 2.3 and 2.4, we need of some technical results, see 
the Propositions 3.1 and 3.2, in section 3, which state the same thesis but use 
some structural and regularity assumptions on the function A. 
Then, these conditions on A, can be removed by a particular compatibility be- 
tween the second order linear term and the first order term in the equation (1.1). 
The proof of these facts is based on the classical property about the preservation 
of the notion of viscosity solution through a global, increasing, change of the 
variable u. 
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3 Proof of the Results 



This part is devoted to the presentation of the technical results which are useful 
for proving Theorems 2.2, 2.3 and 2.4, illustrated in the previous section. 

Proposition 3.1. Consider a viscosity solution u of problem (1.1) valued in a 
bounded closed subinterval I of the domain (a, b), such that u(, i) G W 1,00 (R A ') ; 
with a norm uniformly bounded in the time t G [0,T). Assume a semiconvex 
initial datum uq (see definition 2.1) with a constant Lq > 0. Suppose that, 
H{-,t), w(-,t) are W 2 '°°(R JV ) functions and, f(-,t,-) € W 2 '°°{R N xI), uniformly 
in time. Moreover assume: 

i) A, rje C((a,b))nC 2 (I). 

ii) A < 0, A' > 0, AA" - 2(A') 2 > 7 over I. 

Hi) For every (x,t) G R N x [0,T), w(x,t) G Im{a T (t)). 

Then, there are positive constants C, which depends on sup te j 0T ) || || w 1 . 00 
and Co, which depends on L and Lip(u ), such that 

u(x + h,t) + u(x - h,t) - 2u(x,t) > -C'oe ct \h\ 2 (3.11) 

holds for every (x, h, t) G R N x R N x [0, T). In particular u(-,t) is semiconvex, 
for every t G [0, T). 

The equivalent result for the semiconcave property of the solution, is the 
following. 

Proposition 3.2. Consider a viscosity solution u of problem (1.1) valued in a 
bounded closed subinterval, I of the domain (a, b), such that u(-,t) G W 1, °°(M 7V ) 7 
with a norm uniformly bounded in the time t G [0, T). Assume a semiconcave 
initial datum u (see definition 2.1) with a constant L > 0. Suppose that, 
fi{;t), w(,t) are W 2 '°°(R N ) functions and, f(-,t, •) G W 2 '°°(R N x I), uniformly 
in time. Moreover assume: 

i) A, 77 G C((a,b))nC 2 (I). 

ii) A > 0, A' > 0, AA" - 2(A') 2 > 0, over I. 

Hi) For every (x,t) G WL N x [0,T), w(x,t) G Im{a T (t)). 

Then, there are positive constants C, which depends on sup te [ 0;T ) 
and Co, which depends on L and Lip(u ), such that 

u(x + h,t) + u(x -h,t)- 2u(x,t) < C e ct \h\ 2 (3.12) 

holds for every (x, h, t) G M. N x M. N x [0, T). In particular u(, t) is semiconcave, 
for every t G [0, T). 
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Remark 3.3. We can remark the analogies between the two formulations 3.1, 
3.2. Actually the hypothesis A' > 0, is the same in both the Propositions. This 
can be interpreted as a consequence of the comparison principle, which requires 
a monotonicity for the hamiltonian with respect to u, (see [14], pag. 462). While 
of course, the others expressions in ii) have exactly an opposite sign, showing a 
perfect symmetry between the two formulations. 

The proof of Proposition 3.2 follows by argues as Proposition 3.1 and then 
it is omitted. 

The proof is based on the same technique proposed in [21], which uses a par- 
ticular test function. Nevertheless the nonlinear part in the equation (1.1) does 
not allow to close the proof as in that work. Hence, we turn to a semiconcavity 
property of the nonlinear part with respect to (u, Vu), and to a monotonicity 
with respect to the unknown. 

Remark 3.4. We observe that if a function g is a semiconvex as in the Defini- 
tion 2.1 and is a Lipschitz continuous function, then, for all x, y, z, the following 
inequality holds, 

g(x) + g(y) - 2g(z) > -L(\x - z\ 2 + \ y - z\ 2 ) - Li P ( g )\x + y- 24 

Where Lip(g) is the Lipschitz constant of g. 

Proof of Proposition 3.1. Consider the function v = ue~ ct , where C is a 
nonnegative constant. Then v is an x-Lipschitz continuous function uniformly 
with respect to the time t G [0, T). Moreover v is a continuous viscosity solution 
of the following equation 

d t v - hr{aa T V 2 v) + {ft, Vv) + \{ve Ct )e Ct \a T Vv\ 2 

+r){ve ct )(cr T v,w) + e - ct f(x,t,ve ct ) + Cv = 0, (3.13) 

where (x, t) G x (0, T), and v = u - We now will go to prove that, under 
the assumptions made on A and 77, we have the semiconvexity of the function v, 
and therefore the semiconvexity of the function u. In particular we are going to 
show that, 

v{x,t)+v(y,t)-2v(z,t) > -M(\x-z\ 4 + \y - z\ 4 + \x + y - 2z\ 2 )? (3.14) 

for every x,y,z G R N , t G [0,T), where M = & max(X , Lip(u ))- By 
the assumptions on u, the initial datum is Lipschitz continuous, therefore this 
constant is well defined; this obviously yields the assertion on v, by pluging 
x = z + h and y = z — h. It is easy to see that the above inequality is equivalent 
to the following one: 

v(x,t) +v(y,t) - 2v(z,t) > 

-M[S + z\ A + \y - z\ A + \x + y - 2z\ 2 )] , (3.15) 
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for every 5 > 0, V x,y,z € M. N . 

Hence, fix e, S, 7 > 0, and consider the following test function 

*0, y, z, t) = v(x, t) + v(y, t) - 2v(z, t) + M [S + ^(\x - z\ A 

+ \y -z\ 4 + \x + y- 2z\ 2 )} + e\x\ 2 + e\y\ 2 + e\y\ 2 + (3.16) 

defined for (x,y,z,t) e U = R 3N x [0,T). The assertion (3.15) is equivalent 
to proving that, for every 6, 7 > 0, there exists e = e (5,j) > 0, such that for 
every < e < eo, the following holds: 



inf*>0. (3.17) 

u 

Actually if (3.17) holds, then fixing a point over U, we can send e to zero in 
the inequality 'f > 0, and then sending also 7 to zero we obtain (3.15). Thus 
we limit us to consider the assertion (3.17). 

We assume as usual that (3.17) is false, and will get a contradiction. Therefore 
exist So, 7o > and a sequence Sj — * 0, as j — > 00, such that 



inf*<0, (3.18) 

u 

with S = So, 7 = 70, and e — Ej, for every integer j > 0. If we consider a 
minimizing sequence {xk,yk,Zk,tk) € U f° r ^- By (3.18), the definition (3.16) 
and the boundness of v, we see that (xk,yk,Zk) must be bounded, so we can 
extract a convergent subsequence, which converges to some point (x, y, z, t) € U, 
which, by the continuity of v, is a global minimum point for ^> over U ; moreover 
t (which is obviously less than T) is strictly 

positive. In fact if t = 0, by (3.18) and the Remark 3.4, we obtain, 



> V(x,y,z,0)> 
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> uo(x)+u (y)-2u Q (z) + M[5+j:(\x-z\ 4 + \y-z\ 4 + \x + y-2z\ 2 )] > 

> -L (\x - z\ 2 + \y- z\ 2 ) - Lip(u )\x + y- 2z\ + 
+M[S +](\x- zf + \y - z\ 4 + \x + y - 2z\ 2 )] > 

> {-L a + ^=M}(\x-z\ 2 + \y-z\ 2 ) 

+ [-Lip{uo) + ^=M]\x + y-2z\ > 0. 
V3 

where in the last two inequalities we have used the following relations 



1 



'xl+xl+xj > — =(xi + x 2 + xa), Vii,i2, i3>0, 



x 1 +x 2 > 2xiX2, Vxi,.T2G 
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Hence by the previous contradiction we deduce that t > 0. So the minimum 
point is an interior stationary point of \&. Setting the function 

g(x,y,z,t) = —v(x,t) — v(x,t)+2v(z,t), (3.19) 
*(x,y,z,t) = M[S+^(\x-z\ 4 + \y-z\ 4 + \x + y-2z\ 2 )] 

+s\x\ 2 + e\y\ 2 + e\y\ 2 + ^ (3.20) 

we have that g — $ = — "J has a global interior maximum point at (£,t) = 
(x, y,z,t), therefore we can apply the classical Theorem about the maximum 
principle for semicontinuous functions of M.G. Crandall and H. Ishii, in [9], to 
deduce that for k — | > there exist (bi,Xi) e R x S N , for « = 1,2,3, such 
that, 



(-&;, e P 2 ' Ui, i = l,2 (3.21) 
(|,|,^ 3 )eP 2 %, (3.22) 



and, if denotes the null N x N matrix, we have 



Xi \ 

X 2 < V 2 $(£,t) + K [V 2 $(£,t)] 2 , (3.23) 

o o x 3 J 

b 1 +b 2 + h = (3.24) 

Where for simplifying notations we have set Vi,V2,Vs for v(x,t), v(y,t), 
v(z,t), respectively, and in a same way <!>!, $ 2 , for the partial derivatives of 
$ evaluated at the considered maximum point. Now we compute the derivatives 
of $. Set 

P = ^f(x-z)\x-z\ 2 (3.25) 


q = 2 -^L(y-z)\y~z\ 2 (3.26) 

m = ^f(x + y-2z). (3.27) 


Then 



9M0 = p^>£ 

$1 = 2ex + 2p + m 
$ 2 = 2ey + 2q + m 

$3 = 2ez - 2p - 2q - 2m. (3.28) 
Moreover, if / denotes the N x N identity matrix, we have 
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(3.29) 



By (3.21), and the equation (3.13), we have 



h < -tr(aa T X,)-(^,^)+AKe c? )e c V T ^| 2 

-r 1 (v l e c ' 1£ )(a T <i> I ,w I ) + Cv t + f %1 i = 1,2 

b 3 < ^tr(aa T X 3 ) - ( M3 , $3) - 2A( U3 e c? )e c? |V$ 3 | 2 
- V (v 3 e Ci )(a T ^ 3 , w 3 ) - 2e- c? / 3 - 2Cv 3 . 



(3.30) 



Where we have used the same notation of (3.21) for denoting the functions 
fi, w, /, and we have also omitted the dipendency of the matrix a by t. Adding 
inequalities (3.30) and using (3.23), we compute 



dMlt) < -tr(aa T {X 1 +X 2 + X 3 ))-[( t i 1 ,$ 1 ) + ( f i 2 ,<S> 2 ) + (fx 3 ,<S> 3 )] + 



+ [A(« 1 e c ' ? )e c 'V T *i| 2 - *)(fie Ci )(ff T $i, ttfi) 



+\(v 2 e ut )e ut \a^ 2 \ A - V {v 2 e c '){a T ^ 2 ,w 2 ) 
-2A(« 3 e c? )e c? |ia T $ 3 | 2 - ^e^)^^,^)] + 
+£7(«i + v 2 - 2v 3 ) + e- c \h + f 2 - 2/3]. 



(3.31) 



Now we use (3.28) and (3.29) for estimating the single part in the brackets 
of the inequality (3.31). Let be the I column of cr, and define 



T {i) 

r(0 



for Z = 1, . . . , d, then by (3.23), we have 
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tr(aa T (X 1 +X 2 + X 3 )) = ^ a®) + {X 2 a^ l \ a®) + (X 3 a® , <J (i) » 

1=1 

d 

< ^((V 2 $(C,i)sW,EW) + 
1=1 

K|V 2 $(£t)£ W | 2 ). (3.32) 
Using (3.29), it is easy to compute that 

V 2 $(£t)£ W = 2e£<'>. (3.33) 
So introducing (3.33) in (3.32), we deduce 

d 

^(^(Xi+Xa + Xs)) < 18e^|cr (/) | 2 = O(e). (3.34) 

i=i 

By condition (3.18), we have 



< O(Ve), 

^[(^(M 1 + M 1 ) + ^H 2 ] + M<5 < 2v 3 - Vl -v 2 . 



(3.35) 



Therefore, using the notations introduced in (3.25)-(3.28), (3.35) and the 
regularity assumptions on /i, we obtain 



+ (A*2,$2> + (M3,*3>l < 

< \(m,2p + m) + (fi 2 ,2q + m) - (fj, 3 , 2p + 2q + 2m) \ + 

+2e\{fj,!,x) + (^2, V) + (M3,2)| < 

< |2(/ii - /u 3 ,p) + 2(^2 - + (A»i +M2 - 2/x 3 ,to)| + O(Ve) < 

< sup ||yu(t)|| w2 ,oo [2(—) 1 3(\p\i + \q\i) + 
te[o,T) ZM 

+O l (bl l + M 1 ) + ^H¥(^H)] + o(Vi) < 

1 s s 

< t sup [4(2^3 - «i - «2) + ^[(^^(H 1 + kl 1 ) + ^l TO l 2 ] + 

+^H 2 ]+o(Vi)< 

< 6 sup \\n(t)\\ W 2,oo(2v 3 -v 1 -v 2 )+0(s/I), (3.36) 

t€[0,T) 

where in the last inequality we have again used (3.35). 
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For estimating the nonlinear part we consider the function 



G(u, 9) = \(u)e ct \<r T 9\ 2 - t 1 {u){cj~ t 9, w 3 ), (3.37) 
which depends on («, 9) G [a, b] x R w , and for every u\, u 2 , 9\, 62, define 

A 2 G(u u u 2 , 9 U 6 2 ) = G(«i, 0i) + G(«2, fe) - 2G(^p 2 -, (3.38) 

Let ?7j denotes the value of rj at the point Vie c *, for i = 1,2,3, then using 
notations (3.37), (3.38) and (3.25)-(3.28), we have 

A(w ie c *)e c V T $i| 2 - n(v 1 e ct )(a T 9 1 , wi) + \{v 2 e cft )e cft \(j T ^ 2 \ 2 

-v(v2e cl )(a T <£ 2 ,w 2 ) - 2A(« 3 e C? )e c? |ia T $ 3 | 2 - r 1 (v 3 e cl )(a T <f 3 ,w 3 )] = 

= A 2 G(v ie c \ v 2 e c \ 2p + m,2q + m) + 2 [G ^ 1 + ,P + q + m) - 

-G{v 3 e ct ,p + q + m)] + [ m (a T (2p + m), w 3 - 101) + 

+ m (a T (2q + m),w 3 - w 2 )] +0(Ve). (3.39) 

In the last passages we have used the inequalities (3.35) for estimating the 
residual terms which involve e. 

r]i(a T (2p + to), w 3 - toi) + r] 2 {a T (2q + m),w 3 - w 2 ) = 

= [m - ^{a 1 {2p + m),w 3 -wi} + [r] 2 - ri 3 ](<r T (2q + m), w 3 - w 2 ) 
+2% [(o- T P, w 3 - »i) + (o- T <?, w 3 - w 2 )] 

-77 3 ((t t to, wi + w 2 - 2w 3 ). (3.40) 

Hence, using the Lipschitz regularity of the function u, the regularity of w, 
jointly with (3.35), we deduce, 

[Vi -r l3 ]{a T {2p + m),w 3 -w 1 ) < (3.41) 

6 1 4 

< h'||oo|k T ||oo SUp ||w(t)|| W 2,oo [4||u|| 00 (— — )3|p|3 + 

i£[0,T) ZM 

+ sup \\u(-,t)\\ w i,co(^-)i\p\i\m\\ < 
te[o,T) ZM 

< IIVI|oo|k T ||oo SUp Ht)||w».»[4|H|oo(^T)»|p|* + J SUp ||u(t)|| W i,o. 

te[o,T) ZiW z te[o,T) 



((A)i| p |l+ * | TO |2)i < C r 1 [i(A : )i|p|i + I * | m |2)i. (3.42) 

2M J 11 2M 1 1 yj ~ L 2 V 2M ; 11 42M 1 
In the last inequality we have used the notation 



Ci = 2||?/|| 00 ||cr T || 00 sup ||w(f)|| W 2,oo max(8||tt|| 00 , sup ||u(t)|| w i,oo) (3.43) 
te[o,T) te[o,T) 
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We can repeat the argument for estimating [772 — Vs] { a (2<7 + m), w 3 — w 2 ) , 
obtaining, again using the second relation in (3.35), the inequality 



[Vi - %](c T (2p + to), w 3 - wi) + [772 - r] 3 }(a T (2q + m),w 3 - w 2 ) < 

~ Cl ^ 2F )l(b|1 + + ^2F |m|2) ] " ° l{2V3 ~ Vl ~ V2) (3 ' 44) 
By (3.35) and the following inequality, 



lm+B -H<,i W )ik-g)U»)i + «n- 

then it is easy to estimate the last two terms in the brackets [•] in (3.40), 
through the expression C 2 (2v 3 — Vi — v 2 ), where 



C2 = 5||'7l|oo||c rT ||oo sup ||«;(f)|| W 2,oo. (3.45) 
te[o,T) 

Now we proceed with the estimates of the first two terms in (3.39). We 
observe that, by the assumption ii), we can consider the positive constant, 



n'\UM 2 

2 4 min A' 

and also we can write 



1 llr/ll 2 \\w\\ 



d u G(u, 9) > -C 3 e' ct , V («, 6) e [a, b] x R N . (3.47) 
Therefore by (3.47) and 2v 3 — vi — v 2 > 0, the following holds, 



G{ ^ 1+V 2 2)e ,p + q + m)- G(v 3 e ct ,p + q + to) < C 3 (2v 3 - Vl -v 2 ). (3.48) 
By Hi), we can set w 3 — <r T b 3 , for some iV-dimensional vector 63. Moreover 

8 2 U G = e c h"\<j T 6\ 2 - V "(a T e,w 3 ), 

d 2 G = 2e c ~ t \aa T , (3.49) 

where J = 2e c *A'fl - rfb 3 . For every X — (k,h) e R x R N consider the 
orthogonal matrix A such that Aaa T A T is diagonal, with entries Si > 0, i = 
1, . . . , N. Set h = Ah, 7 = AJ, and define 



c,_ v 1 hWm (3 . 50) 



4 min(2^ - A") 2 max(A) 
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by (3.49), ii), and (3.50), we have 



([VG-C<e-<*[ ' ° ))X,X) 





{d 2 u G - C 4 e- C ')k 2 + ^{2e c ' i \h Z l + 2kJ i h i )S i < 



N 



N _ C 1 

< [%G-C 4 e-<*-'E^S i \1 i \*]h> = 
i=i 

_ P ~ct 

= [e ct A> T 0| 2 + r/'^O, w 3 ), -C A e~ ct - —- \2e ct \'<7 T 6 + V 'w 3 \ 2 ] k 2 



2 , f„/t f)\'J\/„T/ 

_ i i ( ) (7 ■ . v . . . . 

Ct 

W\ 2 \w 3 \ 2 -C,e- Ct ]k 2 < 



2X 

4min(2^ - A") 

P -ct 

-2^llVll a oolHl2 o -fte-«]*» = 0. (3.51) 



Set the constant 



where [•] + denote the positive part of a real number; therefore by (3.51), 
(3.25), (3.35), the Lipschtiz regularity of u, and using the Young inequality 
with exponent 2 we can write, 
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A 2 G(v 1 e Ct , v 2 e c ~\ 2p + m,2q + m) 



A 2 (G - C 4 e- C? y )(vie c \ Vl e c \ 2p + m,2q + m) + 

1 /~i Ct\ 2 2|2 ^ 

+ -C 4 e |w 2 -Wil < 

ci\„, „, i2 



l 



= -C 4 e~ ct \u-i - u 2 \ 2 < 
4 

< ^C 4 (\ Ul ~ u 3 \ 2 + \u 2 - u 3 \ 2 ) < 

< ^^11^)1,^ (i[(A)l (N l + + ^)< 

< ^4 1 sup ||u(t)||^,o. ((2«s - «i - «a) + SM[^ - 1]+) < 

< C 4 C 5 {2v 3 -v 1 -v 2 ). (3.53) 

Introducing estimates (3.40), (3.41), (3.44), (3.48), (3.53), in (3.39), we fi- 
nally obtain 

[A(i; ie c? )e c V T *i| 2 ~ f)("ie C ')((T T $i, Wi) + \{v 2 e c ^)e c '\a T ^> 2 \ 2 

-v(v2e C ')(<J T <i> 2 ,w 2 ) -2X(v 3 e c? )e c ^a T ^ 3 \ 2 -7 1 (v 3 e c ')(a T ^> 3 ,w 3 )] < 

< (Ci + C 2 + C 3 + C 4 C 5 ){2v 3 - Vl - v 2 ) + 0(Ve). (3.54) 

Consider now the last term in (3.31), by the regularity assumptions on the 
function /, the same argument of (3.53) used for estimating \u\ — u 2 \ and again 
(3.35), we have, 

/1 + /2-2/3 < t sup ) ||/(i)|| w2 ,oc[[(^)f(b|l + Mt) + (^) 2 H 2 ]U 

+\ui - u 3 \ 2 + \u 2 - u 3 \ 2 + \m + u 2 - 2u 3 \] < 

< sup \\f(t)\\ W 2^{2v 3 -v 1 -v 2 + 5M[-^-l} + 
te[o,t) 4M 

+2 sup ||u(t)||wi,ooC 5 (2u3- vt -v 2 ) + e ct (2v 3 -v 1 - v 2 )} < 
te[o,T) 

< C 6 {2v 3 -v x -v 3 ), (3.55) 
where we consider 

C e = sup ||/(t)||w2,=c (l + [-L- l] + + 2 sup KOIIwLooCs + e^) (3.56) 
te[o,t) \ 4iW te[o,T) y 
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Finally using (3.34), (3.36), (3.54), and (3.55) in (3.31), we have 



< [6 sup \\n(t)\\ W 2,oo+C 1 +C2 + C 3 + C 4 C 5 
1 te(o,T) 

+C 6 e- C * - C](2v 3 - vi - w 3 ) + 0(>/e). (3.57) 

From the definition (3.56), we see that C§e~ ct is bounded as a function of 
C > 0, so if we choose C sufHcently great, we obtain a contradiction letting 
e — > 0. This prove the result. 



Now we use these Propositions to eliminate the conditions on A and to ob- 
tain the assertion of the Theorems 2.2 and 2.3. 

Proof of Theorem 2.2. We go to build a change of variable such that 
the new differential equation will have the required structural properties of 
Proposition 3.1. 

Let c, A be respectively the infimum of the closed interval I where u take its 
values, and the primitive of A with A(c) = 0. Then consider the solution Q = 
Q{t) of the following ordinary Cauchy problem 

^=cxp(4V7TT+2A(Q)), , 

Q(0) = c. { *- 0S> 

The problem admits an increasing local solution. Moreover by the continuity of 
A, Q is also C 2 . We prove that Q maps the interval I. Consider the following 

cases 

Case 1, (b = oo). Let [0, r*) be the maximal interval of existence for Q, and 
denote with Q the limit for r — > oo. If t* = oo, then Q = oo. Actually by 
the equation (3.58), we have, 

Q'(t) > exp ( inf_ A(u) + 4) , V r > 0. (3.59) 

ue[c,Q] 

So integrating (3.59) from to t > and letting r to infinity we obtain the 
assertion. If r* < oo, then for definition of maximal interval Q blow-ups 
at t* . Otherwise since that A is defined in (a,oo), the solution Q could 
be extended. 



Case 2, (b < oo). Consider again the maximal interval of existence. Then if 
t* — oo and Q is strictly less than b, then, again using (3.59), we obtain 
a contradition. If r* < oo, then and Q < b, then the solution Q can be 
extended bacuase A is continuous in (Q,b). 

In each case the function reachs b in the limit sense; in particular Q maps 
I. Moreover Q can be defined in an open interval V C (— £o,oo), £o > 0, and 
/ C Q(V) C (a, b). By the increasing property of Q, it admits a C 2 (Q(V); V) 
inverse, which we donote P. We use Q as a transformation for a global change 



16 



of the variable u . The function r = P o u, is a bounded, i-uniformly Lipschitz 
continuous viscosity solution of 

d t r hriaa^r) + (/x, Vr) - 

+t 1 (Q(t))(<j t Wt,w) +f(x,t,Q(r)) = 0, (i,t) el" x (0,T), (3.60) 

where the initial datum is To — P o u , and which takes values in the closed 
interval P(I). It is easy to verify the structural hypothesis i), ii) of Proposition 
3.1, where A and r\ are subsituted, respectively, by the functions — (1 + r)~2 and 
i] o Q 1 over the interval V, with regularity properties over P(I). So applying 
Proposition 3.1, we deduce that there exist positive constants C,Kq > 0, such 
that, 

t(x + h,t) + t(x -h,t)- 2t(x, t) > -e ct K \h\ 2 , V x, h e R N , t G [0, i){3.61) 

where the constant Kq depends on P and, on Lq and Lip{u^). Therefore for 
every x, h <G R N , t e [0, t), and some s + , s~ e [0, 1], denoting 

t+ = r{x,t) + s + (t(x + h,t) - r(x,t)), 
t~ = r(x,t) + s~ (t(x — h,t) — r(x,t)), 
t = r(x,t), (3.62) 

we can write, 

u(x + h,t) + u(x -h,t)- 2u(x,t) = Q'(t + )(t(x + h,t)- r(x,t)) 
+Q'(t-)(t(x - h,t) - T(x,t)) 
= Q'(t){t(x + h,t) + t(x -h,t)- 2t(x, t)) + s+Q"{t ++ ){t{x + h, t) 
-T 2 (x,t) + s~Q"(t )(t(x -h,t)- r 2 (x,t) 

>-[K Q C o + Co}\h\ 2 . (3.63) 

for some t ++ S [min(r, r + ), max(r, r + )], r G [min(r, r _ ), max(r, r~)], where 
Coo is a positive constant depending on Q', while Co depends on the Lipschitz 
constant of the solution u and on Q" . This prove the assertion of the Theorem. 



The equivalent result for the semiconcavity property can be then obtained 
with same arguments. So we limit us to give some outlines in the following proof. 

Proof of Theorem 2.3. By the same notations used for proving Theorem 
2.2, we observe that, choosing the increasing transformation u = Q(t), where 
Q, is implicitly defined as the solution of the oridnary Cauchy problem, 

J ^-cxp(-^ T (r + l)'+ 1 + 2A(Q)), 

I Q(0) - c, > 

where /, is choosen as bigger than 3, then we obtain, as in the semiconvexity 
case a new equation for r, which satisfies the structural hipothesis required for 
applying Proposition 3.2. 
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As an immediate consequence of Theorems 2.2 and 2.3, and by the definition 
of the space W 2 '°°, it follows the proof of the Theorem 2.4. This conclusive fact 
allows us to have a second-order regularity result for the solution of problem 
(1.1), with regular initial data. 



4 Time regularity 

In this section we use the spatial regularity of the solution u, obtained through 
the Theorem 2.4, to prove the time regularity of it. 

The result which we show here, see Theorem 4.1, is not stated for second order 
Hamilton- Jacoby equations which have the structure of (1.1), since the lack 
of regularity in the spatial variable. So we present it as a possible interesting 
extension of previous works, in the framework of the viscosity theory. 
In order to simplify the notations, in the following, for every function g defined 
in x [0, T) x (a, 6), and T > h > 0, we set gh(x, t, u) = g(x, t + h,u), where 
< t < T — h, and x £ M. N , u e (a, b). Moreover, if g is a Lipschitz continuous 
function over [0, T), uniformly with respect to the other variables, then we shall 
denote L(g) the constant defined as: 

L(g)= sup \g(x,t,u)-g(x,s,u)\ 



t,*e[o,T) ^~ s| 

t ^ s 
(x,u) eU N x (a,b) 

Theorem 4.1. Let u{t) e W 2 '°°(M. N ), for every t € [0,T). Assume that 
a, /i, w, f are Lipschitz continuous functions of the time, uniformly with re- 
spect to the other variables, and ji, w, f are bounded, then u(x, ■), is a Lipschitz 
continuous function uniformly in R N . 

This result is based on the next Lemma, which shows how the solution u 
tends to the initial datum, when we send the time at zero. 

Lemma 4.2. If m € W 2,oc 7 and fi, w, f are bounded, then there exists a 
constant Co > 0, such that 

\u(x,t) - u (x)\ < C t, (4.65) 

holds for every (x, t) G M. N x [0, T). 

Proof. We limit us to prove the inequality, 

u - M > -C t, (4.66) 

choosing, 

C = snp{H(x,t,u,p,X) : (x, t, u) e M. N x [0, T) x I, \p\ < ||Vuo||oo, 

11*11 < l|V 2 u ||oo}. (4.67) 
The other follows with same arguments. Consider the function, 

$ 7 , £ (x, t) = u(x, t) - u £ (x) + C t + e\x\ 2 + (4.68) 
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For every 7, e > 0, [x, t) e WL N x [0, T). Where ■Ug is the convolution between the 
standard mollificr and u . The inequality (4.66) is satisfied if for every 7 > 0, 
there is e > 0, such that, for every < e < £0, 

inf $ 7 s > 0. (4.69) 

l»x[0,T) 

Actually if (4.69) holds, setting a point (y,s) el ff x (0,T), we have 

u{y, s) - u e {y) + Cos + s\y\ 2 + > 0. 

So we send 7, e to zero and obtain (4.66). 

Suppose that (4.69) would be false. Hence by the same observations used for 
proving the Proposition 3.1, there is 70 > 0, and a sequence £^ — > 0, as j — > 00, 
such that <J> 7j£ has a global minimum point, (x, t), and <f> 7i£ (x, t) < 0, for 7 = 70, 
e = Sj. If i = 0, then 



7 

> $ 7 , e (ar, 0) > -£||u || w i,~ + ^- 



Hence for large j, i > 0, and the minimum point is an interior point; moreover, 
u o~ C t — s\x\ 2 — -j?— £ is a test function for V 2 ~u at (x,t). By (4.68), we deduce 
that 

e\x\ = y/e 2 \x\ 2 = 0(ei). (4.70) 

Moreover by the regularity assumption on uq, we have || Vug II 00 < ||V«o||oo an d 
||V 2 it[j||oo < ||V 2 uo||oo, where Vuo, V 2 i*o, denote the weakly derivatives of uq. 
Introducing the derivatives of the test function in the equation (1.1) and using 
the assumption of boundness of the coefficients, we can write 

Co + ^ < H(xXu(x,t),Wu e a (x,t),W 2 u E (x,?)) + O(e^). (4.71) 

By (4.67), letting j — > 00 in (4.71), we obtain the contradiction. 



Proof of Theorem 4.1. Consider u h , for T > h > 0. u h {t) e W 2 ' 00 ^) 
uniformly in time t <T — h, and is a viscosity solution of the problem 

StUh - ^tr(aa T (t + h)\7 2 u h ) + {fi h ,Wu h ) + X(u h )\a T (t + K)Vu h \ 2 
+T 1 {u h )(a T (t + h)Vu h ,w h ) + f h (x,t,u h ) = 0, (4.72) 
in R N x (0, T - h). For every (x, t) G l w x [0, T), define 

u h (ar,t) = u(x,t) + (C + a(t))h, 

a(t) = eBlt ~ l {CoB l + B 2 ). (4.73) 
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Where Co has been defined in Lemma 4.2, while 

Bi = l|A||oo|k T ||L sup IKOIIwloo 
te[o,T) 

+ hl|co|k T ||co|M|co SUp ||«(t)|| w i.co+L(/). 

te[o,T) 

(4.74) 

B 2 = L(aa T ) sup (^ 2 ||u(t)|| W 2.~ + ||A||oo sup || ||^ ll0 o ) 
te[o,T) 1 te[o,T) 

+ sup ||u(t)|| w i,oc (L((T T )||w;|| 00 + £(u;)||cr T ||oo + L(fj,)). 
te[o,T) 

Then u h is a viscosity solution of the following problem 

d t u h - a'h - ^tr(aa T V 2 u h ) + (/x, Vu h ) + \(u h - C h - ah)\a T Vu h \ 2 

+r 1 (u h - C h - aft)(er T W\ w) + f(x, t, u h - C Q h - ah) = 0. (4.75) 

Recall the Hamiltonian H given by (1.2), then 

we observe that by the regularity assumption on u(t) for every time t, we 
have, for every element (s,p,X) G V 2 - ± (x , h), for (x ,h) £ R N x (0,T), 
\p\ < ||«(to)||w 1 .«= 5 ||-X"|| < C|l u (*o)|| w 2 .°° , for some constant C > 0. 

Therefore for every (ar, t) G 1 N x (0,T - ft), if (s,f>,X) G V 2 '+u h (x,t), 
then (s,p,X) G V 2 ' + u(x, t + ft), and |p| < sup te [ T ) ||u(t)|| w i,oo while |X| < 
C sup te [ T \ ||u(t)||w2,°o. The same kind of observation holds for u . 
Moreover using the last observations, and introducing the constants (4.74), we 
can estimate (4.72) and (4.75), obtaining 

d t u h + H(x, t, u h ,Wu h , V 2 u h ) - B 2 h < 0, 

d t u h + H(x, t, u h ,Wu h , VV) + Bt(a + C )h - a'h > 0, (4.76) 

on the domain R N x (0,T — ft), in a viscosity sense. By the definition (4.73), 
we have 

< d t u h + H(x,t,u h ,Vu h ,V 2 u h ) + B 1 (a + C )h-a'h 

= d t u h + H(x,t,u h ,Wu h ,W 2 u h )-B 2 h. (4.77) 

Hence Uh and u h are subsolution and supersolution of the same problem, rep- 
sectively, and by Lemma 4.2, Uh(x, 0) = u(x, ft) < u + Coh = u h (x, 0), for every 
x G M. N . Using the comparison principle, which holds applying the comparison 
principle which we have stated in [32], we deduce that u(x, t+h)—u(x, t) < C'oh+ 
a(t)h, for every (x,t) G M. N x [0,T). Using the same argument, with the same 
function a, and using again Lemma 4.2, we have u(x, t+h) > u(x,t) — Co—a(t)h, 
for each point (x,t) G R N x [0, T — ft). These two inequalities imply that, for 
every x G R N , h, t 2 G [0,T), 

\u{x, s 2 ) - u(x, si)| < C T \s 2 - si|. (4.78) 

Where the constant Ct depends on the final time T, and B\, B 2 Co- 
rn 

By Theorem 2.4 and Theorem 4.1, then it easily follows Theorem 2.5. 
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5 Regularity and the Ito's formula 



In this section, we derive some consequences from the results of the previ- 
ous sections. Let us recall first from [6] a useful result about Sobolev spaces 
W^' p (M iV x (0,T)), i.e.: the space of the functions u, such that u, Vu, V 2 w, d t u 
ei? oc (R"x(0,T)). 

We start by an observation about the Lebesgue points. 

Remark 5.1. If g € Lf oc (M. N+1 ), for every p <G [l,oo), then for a.e. (x, t), 



r^o+ \n r (x,t)\ Jn r (x,t) 

where Q r (x,t) — B r (x) x (t — r 2 ,t), and B r (x) denotes the open ball in with 
ray r, and centre at x. This follows from a version of Vitali's covering Theorem 
by replacing the balls with Q r , see e.g. Remark 1.3.1 in [16]. 

Proposition 5.2. Let p > N > 2, and u E W^ C ' P (R N x (0,T)). Let 

(x,t) be a Lebesgue point (in the sense of (5.1)) of u, Vu (in L 1 ), and V 2 u, 
d t u (in LP). Then as (x,t) — > (x, t), we have 

\u{x, t) — (u(x, t) + d t u(x, t)(t — t) + (Vu(i, i), x — x) 

+ ^(\7 2 u(x,t)(x - x), (x ~ x))) \ < o(\x-x\ 2 + \t-t\). (5.79) 

For the proof we refer the reader to [6]. If the assumptions of Theorem 2.5 
are satisfied, then we can apply Proposition 5.2 to w, obtaining that for a.e. 
(x,t) e R N x (0,T), we have 

d t u(x, t) + H(x, t, u(x, t), Vu(x, t), V 2 u(x, t)) = 0, (5.80) 

where dtu, Vu, V 2 u, represent the wealkly derivatives of u. This regularity of 
the solution u has a useful application when we use the Ito's rule along a stochas- 
tic process. 

As we have noted in the Introduction, the representation of financial derivatives 
can be expressed as a conditional expectation over the probability space of the 
underlying factors that determine the instrument's price. 

In the equality (1.4), we introduce a relation between this conditional expected 
value and the solution of (1.3). This is a consequence of the regularity of 
u = U + h + £, which follows by applying the results of the previous sections to 
(1.3). In the model proposed by X. Gabaix in [13], the measure Q depends in a 
nonlinear way by the price of the Mortgage-Backed Security and its volatility; 
this feature produces the quadratic nonlinear term in the equation (1.3). To rig- 
orously derive the equation (1.4), we need to apply the Ito's formula to u(X t , t), 
where X t is the underlying stochastic process. When u is less regular than C 2,1 , 
a first result was established by Krylov [25]. In a more recent work, Haussmann 
[17] describes a result in this direction, which shows that the Ito's formula holds 
for u G W 2;1,00 (M Ar x (0, T)), provided that it is interpreted appropriately, using 
the generalized Hessian. 

We use the regularity of the solution U, which gives the property (5.80), com- 
bining it with the result of Haussmann to obtain the equality (1.4). 
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In the next we shall assume p € W 2,1 '°°(R Ar x [0, T] and for sake of simplicity 
we limit us to treat in detail the case of a constant matrix a and consider the 
case N > d. 

As we have already pointed out in the Introduction, equation (1.3) is equivalent 
to the general problem (1.1). Actually, after the change u = U + h + £, u solves a 
problem which has the same structure of (1.1), where in particular w — cr T • Vh. 
Hence, following the comparison principle and the Lipschitz regularity for vis- 
cosity solutions of (1.3) proved in [32], and assuming h to be a smooth func- 
tion, as in [32], by Theorem 2.5 is straightforward to assume the existence of 
a unique viscosity solution U € W 2 ^'°°(R N x [0, T)) of the equation (1.3) with 
U(x, 0) = 0, such that U + h + £ > 0. 

Theorem 5.3. Let (ft, P) be a probability space and let (Wt, Ft)te[o,T] be a 
d- dimensional continuous standard Brownian motion over ft and let X t be the 
Ito process defined as the solution of the stochastic differential equation (1.5), 
for some intial datum X £ R N . Define 

a T (T-s)VU(X s ,T-s) 
7s ~ P U(X a , T-s) + hjx7, T-s)+£(T- s) ' a - S - 



W s = W s + 



/ 7(K)dK, 

Jo 



with < p < 1 a parameter. If N = d we do not need any assumption. 
Otherwise for the case N > d, assume that the projection-process n t , defined as 
the projection of X t over the Kernel of cr T , satisfies 

A ^ P (tt s e A\X t ) , A C R m is £ m -absolutely continuous, (5.82) 

for every s > t, where C m denotes the Lebesgue masure over R m and m = 
N - rank(a T ). Then 

, fT 

U{X t ,T-t)=E?[J (t -r(T - s))e-tt r{ - T -^ dK h{X s ,T ~ S )ds\, (5.83) 

holds a.s. for every < t < T. 

We remark that a classical result, the Girsanov's Theorem, state that W s is 
a Brownian motion with the same filtration of W s , see for istance [1] or [23]. 
To prove the relation (5.83) for the solution U to the model (1.3), we shall use 
the results of [17]. Actually the assertion of Haussmann interprets the Ito's rule 
through some processes which substitute the usual derivatives of the function 
evaluated at the process (X t ,t), which are equal if the underlying process be- 
longs to some set of full Lebesgue measure. In order to state our formula we 
need to neglect the term which correponds to an integration over the paths of X 
which fall in a set of null measure. Therefore the assumption (5.82) can be ex- 
plained by that purpose. Actually the degeneration of the volatility matrix plays 
a crucial role: in general along the components of the vector X which correspond 
to the directions of the kernel of cr T , we have almost- deterministic trajectories 
which surely form a set of zero measure; this feature, in general, could not allow 
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to justify the previous assertion (5.83) about U , using the Haussmann's Theo- 
rem, unless the drift coefficient has some compatibility condition with a; in fact 
this feature is expressed through the condition on the projection-process 7r t , see 
the next remark 5.6. 

As we have pointed out in the Introduction, there are some works such as [5], 
where the existence of densities for the solutions of stochastic differential equa- 
tions under the Hormander's condition is proved. In particular Theorem 3.1 
of chapter 3 by [5] gives a sufficient condition, which involves the use of the 
Lie bracket between the set of vector fields represented by the drift coefficient 
and the columns of the volatility matrix, to prove the existence of an absolutely 
continuous distribution with respect to the Lebesgue measure. Although that 
condition can be used in practice to investigate differential equations with a 
complicated form, it does not allow to consider a time dependence of the coef- 
ficients, and it requires also more regularity for computing the Lie brackets of 
the vector fields. Moreover, while for financial purposes the coefficients \i and 
a can be taken as constants, the probabilistic Hormander's condition used by 
D.R. Bell in [5] does not hold in this case. Therefore we propose an alternative 
assumption, contained in (5.82), which covers the cases of interest for our appli- 
cations. Actually let us notice that case N = d, the assumption (5.82) can be 
removed via an approximation of the underlying process X. The proof of the 
following technical Lemma motivates of condition (5.82). 

Lemma 5.4. Let Q be a probability measure equivalent to P over Q and let B 
be a set of zero Lebesgue measure over R N x [t, T], for < t < T. Set 

Z t , T = J p s lB(X s ,s)ds + q s l B {X s ,s)dW s , (5.84) 

where p s ,q s are (J r s )t< s <T-adapted bounded processes. Then, under the assump- 
tion (5.82), it holds 

E?[|Z t , r |] =0. 

Remark 5.5. Let us recall the following important fact. For every integrable 
random variable Ft -measurable there holds 



E? [X] = 





' Y dQ~ 
A dP 


E* 


' dQ~ 
L dP 





where E t [-] denote the conditional mean with respect to the measure P. To 
see this relation between the conditional means taken through two equivalent 
probability measure, we refer the reader to [30]. 

Proof of Lemma 5.4. By remark 5.5, the boundness of the processes p s , 
q s and 

i 

T 2 1 5 

[ q s l B (X s ,s)dW s 



E, 



q s l B {X s ,s)dW s 



< E 



E 



\q s \ 2 l B (X s , s)ds 
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where in the last passage we have used the classical Ito isometry, we are reduced 
to prove 



E, 



Ib{X s , s)ds 



= 0. 



(5.85) 



Case N > d. Let {61, ... ,b m } be an orthonormal basis of the kernel of cr T , 
then consider an N x N invertible matrix M, such that 

M T a = bi, i = 1, . . . , m, 

and define the process Y s — MX S . Therefore 

Y s = (vrf ,G S ) := ((M T )^ n s , G s ) , 



and 



s)ds, vrfeM™, 
(i G (Y s ,s)ds + a G -dW s , G,€R N - 



(5.86) 



Moreover 



^(y,s) = {{ii{M~ l y, T- s),h),... , {^{M~ l y, s), b m )), (5.87) 
and ere G .Mjv_ m ,d(K) is a matrix of rank N — m. Let 







M = 



M : 


\ ... 1 / 



B = MB. 



Noting that B has zero Lebesgue measure, (5.85) is equivalent to 



E, 



0. 



(5.88) 



Fix e > and t + > t/(l — e), then consider the approximating process = 
{^,GD, where 



r M,e 



f s(l-e) 



X)d\ 



(5.89) 



Gl = G t + J° {1 £ \ G (Y x ,X)d\ + a G -(W s 



W t ). 



It easy to see that this process converges in L 1 to Y s , uniformly in time, as 
e — > 0. The JF s ( 1 _ £ - ) -conditional density of the random variable G e s is normal 
with mean 

r s(l-e) 

E G :=G t +J^ f, G (Y x ,X)dX 
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and covariance matrix 



T G (s,t) := (s-t)a G al > 0. 

Consider S > 0, then there exists a countable collection of sets Qj C R N x [t + , T] 
which, without loss of generality, can be supposed to be closed TV+l-cubes, such 
that the projection of Qj over the coordinates x m +i, ... , xn lies in a bounded 
set of A 2 (this is not a restriction bacause we could always take the intersection of 
B with an increasing sequence of sets whose projection over these coordinates is 
bounded), the intersection of a couple of cubes is a set of zero Lebesgue measure 
and 



and 



B n (R N x [t + ,T]) c UjQj 



J2£ N+1 (Qi)<s. 



We proceed at first with an estimate of the expected value (5.88) where B is 
substituted by Qj and Y s by F/ . By the distributional property of the approx- 
imating process we have 

/ ' ••• [19,07,*)] <& = 
Jt+ 

= [ T E t [E s{1 _ E) [l Qj (Y s s ,s)]]ds 
Jt+ 

= [ E t [P((Y s s ,s)eQ j \F s(1 _ £) )]ds. (5.90) 
Jt+ 

We shall use a diagonalization procedure to separate the variables in the follow- 
ing integration. Hence consider the orthogonal matrix H, such that H T T G H = 
Diag{\\, . . . , Ajv-m), where Ai, . . . , Ajv-m ar e the eigenvalues of T G . 
If Qi(a), Bi{a) respectively denote the cube of centre a and side of lenght I > 
and the open ball of radius I and centre a in R N_ra , then again by the properties 
of G% , we have 



P(G e 8 eQi(a)\r a{1 - e) ) < 



< 



L 



i(<0 m ^fdet(T G {s,t)) 



1 



— — —N—rri — — — — 

^2vr(s - t) y/\l---\N 

1 



< 



-N—m 



= n 



QivTv^r( a ) \/2tt(s - t) \J\\ - ■ ■ Xn-tu 
n 



exp (-^{Ta\s,t){x - E G ),x- E G )^j 

(-^ew)* 



WN — m + ^ 



< 



(2WN -m) N - m 



WN — m + ai 



y/2Tr(t+-t)Jdet(a G al) 



C(N,m,t,t + ,* G )£ N - m (Qi(a)), 



(5. 
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where a = H T (a — Eg), and <f> denotes the standard Gaussian distribution. 
Writing Q 3 = Q) xQ^x Ij, for some cubes Q) C K m , Q) C A 2 , and Ij C T], 
(5.90) can be continued using the estimate (5.91) in the following way 



E, 



l Q ,{^)P{Gl 6Q 3 2 |f s(M ) 



< C{N,m,t,t + ,a G )C N - m {<£) 



ds < 

1 Q 1 A 7r s ) 



ds. 



(5.92) 



By the condition (5.82), there exists the density d w (x; r, u>, t) > P — a.s. of 
the measure P(tt~ (-)\Xt) with respect to C m which, by the regularity of X t , is 
a continuous function of r e By this property, (5.87) and denoting M 

the restriction of M T over the kernel of er T , we obtain 



E, 



l Q 1 Ms ) 



<MQ) 

Introducing (5.93) in (5.92), we have 

/ T E t [l Qi (Y a e ,s)] ds< 
Jt+ 



= / ^(x; s(l — e), w, f)da;. 



(5.93) 



= c(iY,™,t,t+,a G )/: Ar - m (g 3 2 ) 



/ / dn(x; 



s(l — e), X t ,t)dxds 



= C(N,m,t,t+,a G ) [ h 

JMQ]xQ 2 ,xIi 



> {y)d-n(x; s(l — e), u>, t)dxdyds, (54).4) 



lMQ)xQp 

By the properties of the collection {Qj}j, we can write 



J, .WO; x Q| x I 3 



\detM\C N+1 (UjQj) < \detM\5. 



(5.95) 



For P — a.s., the measure v{-) = v(-;u,t), defined through the density function 
(x,y,s) i ► lA 2 (y)d 7r (x; s(l — e), w) overtax is an absolutely continuous 

measure with respect to N + 1-dimcnsional Lcbesgue measure. Therefore set 
/3 > and choose 8 so that for every set A with £ A ' +1 (A) < 8\detM\, it holds 
K-^) < Hence using (5.95), we deduce 



/ lB(Y s e ,s)ds 
Jt+ 

< v (yJjMQ) x Q) x Ij^j < 13. 



(5.96) 



Then the assertion in the case N > d follows by the arbitrary choice of j3 and 
letting e — > and then t + — > f in (5.96). 
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Case N — d. Introducing the approximation 



X £ — Xt 



n(X Xl T-X)d\ + a £ (W s - W t ), 



(5.97) 



where a £ — ^ a, as e — > and det{a £ ) ^ 0, then we can repeat the same arguments 
used for the previous statement substituting X £ at the process Y £ of the previous 
case. Therefore we obtain 



/ l B (X £ s ,s)d i 
Jt+ 



< C{a £ )5, 



(5.98) 



for arbitrary S > 0. Hence the limits with respect to e and t + is taken. 



Remark 5.6. The condition (5.82) in some sense in optimally. Actually in the 
case N > d, the strong degeneration of the quadratic form <t<t t could make false 
the property (5.84). For istance if N = 2 and d=l, choosing 

er T = (0, 1), M = 0, t = 0, B = {0} x (—oo,0) x [0, T], 

then with Xq = 0, 

X t = (0,W t ) 



and 



E 



/ l B (X„a))dfi 
Jo 



T 

2>° 



which contradicts the assertion of Lemma 5.4. 



Proof of Theorem 5.3. Let E(i) = cr(T - i), A*°(a;,i) = //(a:,T - i), 
/i° = h(x,T-t) and = £(T-t), r°(i) = r(i) for every G K w x (0,T]. 
Then define the function U°(x,t) = U(x,T — t). U° is a viscosity solution of 
the equation 



1 iy T vr/°i 2 

= r°{U + h )-Th° and U°(x,T) = 0. 



(5.99) 



Moreover the equality (5.99) holds for a.e. (a;,t) € x (0,T) with the weakly 
derivatives of U° . Applying Theorem 3.1 pg. 733, in [17], to the function 
(x,s) i — ► e~$t r ^ dK U°(x, s) relatively to the time interval [t,T) there exist 
processes /3q, (3, a, such that 







e-^ r °^ dK U°(X T ,T) = U°(X u t) + J T [0o(u>,8) + (»°(X s ,s),P(cj,s)) 

+ itr(SS T (s)a(w,s))]ds+ f (3 T {lj , s)Z(s) ■ dW s , a.s. (5.100) 
2 Jt 
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Moreover there exists a set A C 



piV 



(t, T] of full Lebesgue measure such that 



the usual derivatives of U° exist in A and 

e" J ' r ° Wi " (d s U (X s (cu), s ) - r (s)U°(X s (u),s)) , 
e - // <- ( K )^vf/°(I s ( w ), s ), (5.101) 



A)(w,«) 

a(w,s) = e-^ r °^ dK V 2 C/ (X s H, S ), 



whenever (X s (a>),s) € A. Without loss of generality we can also assume that 
the equation (5.99) holds over A. The application of (5.99), (5.100) and (5.101) 
yields 

= U^X t ,t) + £ [f3 (s) + ( f i°(X s ,s),(3(s)) + ^tr(XX T (s)a(s) 

-/3 T (s)£( s ) 7s )] 1 A (X S , s)ds + J 1 a (X S7 s)(3 t {lj, «)£(*) • dW s + Z t , T 
= U°(X U t) + ^ e-h s r°Wd* [ r °( s )(u°(X s , s) + h°(X s , a)) - rh°(X s , s) 
-r°( S )U (X S7 s)]d S + J f3 T {cj,s)X( S )l A {X s ,s)-dW s 



a.s. 



(5.102) 



By the boundncss of the processes (3$, f3, a, we recognize that the remaining 
term Z t T has the same structure of (5.84) in Lemma 5.4, therefore, by the 
definition of Q, (5.81), as a measure equivalent to P and Lemma 5.4, Z t .T has 
a null contidional Q-expccted value. 

Using this remark and taking the conditional expected value in both the left 
and right side of the equation (5.102), we obtain 



U°{X u t) =E? 



^ e-K r °^ dK ( T - r°(s))h°(X s ,s)l A (X s ,s)ds 



a.s(5.103) 



By the full measure of A and noting that the functions h and r are defined 
everywhere and are in particular continuous, the assertion is proved. 



Finally we obeserve that the same conclusions can be obtained for a process 
X t such that satisfies 

dX t - n(X u T - t)X t dt + aiT - t)X t dW t , Xi(0) > 0, i = l,... ,N. 

In this case can is possible to treat the problem using the same previous argu- 
ments simply making the change y — log(x) in the equation (1.3). 
Also if the drift coefficient depends in a linear way by x, and the diffusion co- 
efficient is constant with respect to x, then we could repeat the same argument 
used before to state the regularity of u. Actually, in that case, we have a linear 
increasing with respect to the variable x, only in the first order linear term of 
the equation. 
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